Lattice softening effects at the Mott critical point of Cr-doped V2O3 
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We have performed sound velocity measurements in (Vi_i:Cri:)203 in the vicinity of the critical 
point of the first order Mott transition line. The pressure sweeps at constant temperature reveal 
a large dip in the C33 compression modulus, this dip sharpens as the critical point is approached. 
We do not observe signs of criticality on the shear modulus C44 which is consistent with a transition 
governed by a scalar order parameter, in accordance with the DMFT description of the transition. 
However, the amplitude of the effect is an order of magnitude smaller than the one obtained from 
DMFT calculations for a single band Hubbard model. We analyze our results using a simple model 
with the electronic response function obtained from the scaling relations for the conductivity. 



I. INTRODUCTION 

The Mott transition, which is a metal-insulator transi- 
tion (MIT) without lattice symmetry breaking, observed 
in various correlated electron systems, has been one of 
the major challenges in condensed matter theorjffl. The 
theoretical interest in this phenomenon has been boosted 
in the recent years by the success of the Dynamic Mean 
Field Theory (DMFT) approach of the Hubbard modeP. 
One of the theoretical difficulties concerns the choice of 
the order parameter. Indeed, as the transition does not 
break any lattice symmetry the order parameter is not 
clearly identifiable. It was argued that the MIT should 
be analogous to the liquid-gas transition and belongs to 
the Ising universality class^l. Various quantities express- 
ing the " metallicity" of the system have been employed 
to represent an order parameteJ^''^. 

Another issue is related to the role of the lattice. Trac- 
ing back to the earliest ideas, the electron-lattice coupling 
had often been envisaged as the possible mechanism for a 
first-order MIT. Although the modern DMFT approach 
can account for the behavior of the resistivity in the vicin- 
ity of the critical line within a purely electronic model^, 
lattice degrees of freedom do play a role at the transition 
of real materials, where it is accompanied by a discontin- 
uous volume change^^^^. 

On the experimental side, examples of widely studied 
materials are the Cr doped vanadium oxidt"^ and organic 
charge transfer salts k-BEDT2X^^ In the latter, ev- 
idences for a lattice softening near the critical endpoint 
of the MIT line were found in ultrasound experimental^ 
and by observation of anomalies in the thermal expansion 
coefficientJi^. 

The Cr-doped vanadium oxide (Vi_2;Cr2;)203 is the 
archetype compound where the MIT can be induced by 
doping or applying pressure. Here the transition line 
ends with a second order critical point near to T = 450 
K?^. The criticality at this point was studied by electri- 
cal transport which showed scaling properties expected 
by DMFTGI. However, experimental facts concerning the 



role of the lattice in this transition are rather scarce. 
The information about the lattice response should be 
important for testing and refining realistic models which 
take into account the orbital degrees of freedom. Indeed, 
recent simulations using state-of-the-art LDA-I-DMFT 
techniques indicate that in V2O3, the interplay between 
the correlations and the orbital polarization effects mod- 
ifies substantially the nature of the MIT compared to 
that of a genuine Mott transition in the one-band Hub- 
bard modep3. 

In this work we try to shed a new light on the ques- 
tion of the lattice contribution to the MIT through ultra- 
sound measurements in the vicinity of the critical point 
of (Vi_a;Cr2;)203. We show that this experimental ap- 
proach can address both the questions of the symmetry 
of the order parameter and its coupling to the lattice 
degrees of freedom. 

II. EXPERIMENTS 

All measurements were performed on single crystals 
with a nominal dopingof x = 1.1%, prepared using the 
skull melter technique^^. This Cr concentration ensures 
that the sample is on the insulating side of the transition 
at ambient pressure, but a moderate pressure of a few 
kilobars drives the system into the metallic state. This 
doping also corresponds to the one used for transport 
studies in^. 

In this study we measured the propagation of acous- 
tic waves along the hexagonal c-axis, for both longitu- 
dinal (compression) and transverse (shear) mode. The 
measured sound velocities are proportional to the square 
root of the corresponding elastic constants; for the lon- 
gitudinal mode it is the compression modulus C33, and 
the transverse waves involve the shear modulus C44 (the 
two transverse modes are degenerate for the propaga- 
tion along the c-axis due to the trigonal symmetry of 
V203)^^. As it will become clear later, the choice of a 
high symmetry axis is crucial for the interpretation of 
our data. Here we can measure pure compression and 
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shear modes, whereas for almost any other direction of 
propagation the excited modes would involve a mixture 
of longitudinal and transverse components of the strain 
tensoiESl. 

The crystals were oriented by X-rays, and cut along 
planes perpendicular to the c-axis to an overall thickness 
of approx. 2-3 mm. They were then polished on a lapping 
machine in order to achieve parallel planes on both sides 
of the sample and a surface roughness of less than 1 /im 
(roughly a tenth of the wavelengths used). The velocity 
measurements were performed using the standard tech- 
nique where ultrasonic waves are generated by applying 
short radiofrequency pulses to a piezoelectric transducer 
and a phase sensitive detection is used to measure the 
echo signal^. We note that no frequency dependence of 
the sound velocity was observed within the used range of 
200-500 MHz. 

One of the technical difficulties encountered is that we 
cannot use commercial transducers for sound generation 
and detection, because in high temperature and pressure 
condition transducer bonding proves difhcult to achieve 
and unreliable. Instead, the transducers were directly 
grown on the sample by sputtering: a piezoelectric ZnO 
transducer was grown on top of a thin chromium bonding 
layer and a gold contact layer. The thickness of the ZnO 
film was about 10 /im corresponding to an acoustic reso- 
nance frequency in the 400 MHz range. Since ZnO crys- 
tals grow with their main piezoelectric axis parallel to the 
c-axis, in principle only the longitudinal mode should be 
excited. Though, due to imperfections in the ZnO layer it 
was possible to excite and observe small transverse mode 
echoes on some samples. The two modes having very dif- 
ferent velocities it is in principle easy to unambiguously 
separate them when both are excited (thanks to a very 
broad frequency response of the transducer we could do 
it using very short pulses, less than 100 ns). Very unfor- 
tunately, for T > Tc the two velocities are separated by a 
factor very close to 2 so that the echos of the two modes 
strongly overlap and it becomes difficult to resolve them 
even with the shortest available pulses. Because of this 
accidental condition we could not obtain the complete 
set of data for the transverse mode. 

The high pressure conditions were obtained using 
isopentane as pressure liquid and an externally controlled 
pressure system where the cell is connected by a capil- 
lary tubing to a pressure generator. High temperatures 
are obtained with a small heater inside the cell. Ex- 
periments were performed by sweeping continuously the 
pressure (at a rate of ~20 bar/min) at constant regu- 
lated temperature. Our measurements are limited to a 
pressure of roughly 5 kbar for temperatures in the range 
of T 500 K. These limitations are mainly due to the 
materials used in the pressure plug and electrical connec- 
tions. 



III. RESULTS AND DISCUSSION 

The evolution of the velocity of the longitudinal sound 
waves with varying pressure at several fixed tempera- 
tures is shown in Fig. [T] The data plotted here have 
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FIG. 1: Pressure dependence of the relative velocity variation 
Au/w for the longitudinal mode at various temperatures about 
the critical point. 

been corrected for the usual linear temperature depen- 
dence due to the thermal expansion of the lattic«P3 
Au/w = (u(p, T) - vun{T))/viin{T) where vu„{T) is a 
linear function of T representing the temperature depen- 
dence of the sound velocity observed far from the crit- 
ical regime (here for p ^ 500 bar in the interval 250- 
400 K). Such a linear temperature dependence has been 
observed for various elastic con stants and it is usually 
attributed to anharmonic effect^^SEI] Therefore Fig. [l] 
represents solely the critical contribution to the sound 
velocity. Note that for T ^ Tc the critical regime is span- 
ning a pressure domain large enough to affect the sound 
velocity even at ambient pressure, as can be seen from 
the small negative offset of the curves for temperatures 
close to Tc in figure [l] At temperatures where a scan in 
the full pressure range was taken, the relative sound ve- 
locity shows a pronounced minimum around Pc 4000 
bar. The largest effect is observed around 458 K which 
is close to the critical temperature of 457.5 K estimated 
from transport measurement^. 

The most challenging aspect of these experiments is 
related to the extreme sensitivity of sound propagation 
to the sample integrity (we emphasize that our samples 
are many times bigger than those used for transport mea- 
surements in'^). The sharp volume change at the transi- 
tion (this change amounts to 1% at room temperatur^^ 
most often gives rise to microfractures resulting probably 
from releasing some internal stresses frozen during crystal 
growth. Such fractures spoil or even suppress the acous- 
tic echo signal and eventually make the sample unusable. 
In practice this means that it is forbidden to cross the 
first order transition line during pressure or temperature 
sweeps. For this reason, below the critical temperature 



3 



pressure scans were only performed in the insulating low 
pressure region. Note also that very close to Tc (below 
470 K), the sharpness of the minimum of Av/v depends 
on the sample (c.f. two curves for 458 K), proving that 
it is also very sensitive to the crystal quality and most 
likely the differences in the doping inhomogeneity. In 
fact, in samples showing the sharpest curves the micro- 
fractures often appear even slightly above Tc- Therefore 
in order keep the sample unaltered during the set of pres- 
sure sweeps and to get reproducible results we are led to 
make a compromise and choose a less homogeneous sam- 
ple which shows a more smooth pressure variation near 
Tc- All data shown in figure[l]were obtained on the same 
crystal, except one of the curves at 458 K which comes 
from a sample showing the most pronounced dip. 

We now turn to the analysis of the transverse mode 
velocity, shown in Figj2] The vertical offsets of different 
curves were corrected for the temperature dependence 
by the same procedure as described for the longitudi- 
nal mode. As explained in the previous section, these 
measurements could only be carried for T < Tc where 
the transverse echoes are not hidden by the much larger 
longitudinal echoes, in this region complete scans were 
not performed to avoid the expected harm to the sample 
when crossing the first order transition line. Nonethe- 
less the data taken up to 448 K show the general trend 
which is much different that the one observed for the lon- 
gitudinal mode. In the region closer to the critical point 
all curves show a similar slow upward shift, moreover all 
curves between room temperature and Tc stay parallel in 
the region p ^ 2-3 kbar, with no signs of critical behav- 
ior at Tc- The small increase under pressure is similar 
to that observed at room temperature in pure V2O3, a 
compound which does not exhibit the MIT^H. The over- 
all behavior of the transverse mode is in contrast to the 
data of Figjl] where in the same temperature interval a 
dramatic drop, developing already below 2 kbar, is sharp- 
ening as we approach Tc- While the missing data at the 
critical point does not allow us to consider this as a rig- 
orous demonstration, we think that this finding gives a 
reasonable evidence that the transverse mode does not 
exhibit any critical softening close to the metal insulator 
transition. 

For an Ising-class MIT, the absence of criticality of the 
shear modulus can be inferred from simple symmetry ar- 
guments. A quite general approach to the effects of a 
phase transition on the elastic constants was developed 
by RehwalcP^ in the framework of the Landau theory, es- 
pecially to study soft modes of lattice vibrations observed 
in structural transitions. In this approach the Landau 
free energy functional contains coupling terms which rep- 
resent the interaction between the order parameter 77 and 
the strain tensor e. The general expression can be written 
as a linear combination of products of increasing powers 
of some components of both the order parameter 77 and 
the strain tensor e: 
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FIG. 2: Pressure dependence of Av/v of the shear mode at 
various temperatures 



where the nonzero coefficients a, /3, 7,... are selected by 
symmetry. 

It has been shown that a singular softening of an elastic 
constant such as the one in Fig. [T] requires the existence 
of linear coupling between 77 and e. The possible combi- 
nations can be enumerated by symmetry considerations. 
In the language of group theory only the components of 
1] and e belonging to the same irreducible representation 
of the point group symmetry of the disordered phase can 
be coupled. For a scalar order parameter this means 
that allowed combinations of the strain tensor elements 
should form the basis functions of the one dimensional 
representation Aig. In our case, this requirement is ful- 
filled for the ezz component involved in the compression 
mode, since it is invariant under the operations of the 
trigonal symmetry group. On the other hand one cannot 
construct an invariant from a linear combination of ezx 
and ezy, and therefore there is no linear coupling allowed 
for the transverse mode. Our failure to detect a critical 
contribution to the shear mode velocity is therefore con- 
sistent with the fact that the Mott transition is expected 
to be governed by a scalar order parameter. 

In the following we will try to quantitatively analyze 
the data of Fig{l] First note that the maximal ob- 
served amplitude of the decrease of the relative speed 
of sound Av/v ~ 2% is one order of magnitude lower 
that the value reported for organic compounds. It is also 
smaller by the same factor when compared to the theo- 
retical values given by Hassan et al^ based on DMFT 
calculations for a single band Hubbard model (e.g. at 
T = 1.27Tc w 580K these calculations give a 3% effect). 
On the other hand there is a qualitative agreement on the 
pressure dependence, namely the width of the dip which 
is of order of 2 kbar for T - 500 K. 

Here we will reconsider the same model in a slightly 
different way. The compressible Hubbard model is based 
on the work of Majumdar et al^, who presented the 
first approach to include lattice effects in this theoret- 
ical framework and treated it in the simplest approxima- 
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tion where all phonon excitations are neglected. They 
take into account the dependence of the free energy on 
the unit cell volume V through the expression F{V) — 
Fe{D{V)) + E{V) where is the electronic contribution 
depending on V through the variation of the bandwidth 
D, and E{V) is the elastic lattice energy. The lattice 
stiffness (compression modulus) 

K = -VdP/dV Vd^F/dV^ 

can then be expressed as a sum of lattice and electronic 
contributions: K = Ki -\- Kg,. The latter is negative and 
is usually written as 

K,^-V{dD/dVfxei 

where the electronic response function is defined as 

Xei {T,D{V)) = ~-d^FjdD^ 

The purely electronic Mott transition is marked by a di- 
vergence of Xei at T = Tel, however in the compressible 
model the instability comes from a divergence of the com- 
pressibility 1/K oo. At this point a first order MIT 
accompanied by a discontinuous volume change occurs 
to avoid the unphysical region where < 0. Since this 
happens for a finite value of Xeh the true critical tem- 
perature Tc > Tel- The sound velocity is proportional 
to K-^/^ therefore it was predicted that it would van- 
ish as (T — Tc)^/^^. However, this expansion is derived 
assuming K close to zero, i.e. \Ke\ ^ Ki. Obviously, 
with Av/v of order of 1% we are far from reaching such 
a regime. A more realistic scaling will be obtained as- 
suming \Ke\ « Ki in the domain studied here. In this 
limit the sound speed variation obviously follows that of 
Xei- IiJ^ XeZ was determined from DMFT calculations. 
Here we adopt a different approach using the available 
experimental data. Limelette et al^ argued that in the 
vicinity of the critical point the excess conductivity in 
the metallic state a — acrit can be regarded as an or- 
der parameter. It was shown that this quantity verifies 
a universal scaling relation, which, for T > Tc is writ- 
ten as h^/^f{ht - -f5{5 - 1)) where t = {T - Te)/Te, 
h = {P — Pc)/Pc, and the exponents were shown to be 
close to the mean- field values 7 = 1 and (5 = 3. We 
have calculated Xei = da/dh from the experimental scal- 
ing function / given irP. The apparent symmetry of our 
experimental curves allows us to suppose that the same 
scaling in function of can be used on both pressure 
sides, which is in general expected for T > Tc if is 
considered as a conjugate field (note that the asymmetry 
in the DMFT simulation^ii' is much less pronounced for 
T > Tc). The results of simulations are shown in Fig- 
ure |3] In contrast to the DMFT calculations, we have to 
treat the amplitude of Xei as an adjustable parameter. 
Therefore Tei cannot be independently determined. In 
our simulations we have taken Tei = 450K for which the 
best agreement of the pressure dependence can be ob- 
tained at least above 470 K. Our curves go to zero about 
451 K whereas the value iiP gives Te-Tei « 1.4% w 6K. 



This discrepancy comes again from the overestimation of 
Av/v and Xei by the same factor. Our approach, based 
on the assumption that conductivity is directly related 
to Xeh means that we take into account only the effect 
of the electronic degrees of freedom on the lattice, but 
neglect the feedback of the lattice on the electronic de- 
grees of freedom. This simplification is now justified by 
the very small influence of the lattice on the transition 
temperature. 

To check the effect of the sample homogeneity we also 
calculated the theoretical values of v assuming a gaus- 
sian distribution of Pc (dashed curves in Figure |3| . To 
reach the experimental values of the minimum at 458 K, 
it is necessary to assume a 20% gaussian width, an un- 
likely high value which is clearly incompatible with the 
sharpest curve at this temperature. We may speculate 
that either there is a mechanism precipitating the tran- 
sition already at very low values of \Ke\/Ki, or at some 
point the divergence becomes too narrow to be observed 
in real samples. 
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FIG. 3: Comparison of simulations for Tei ~ 450-R' (solid and 
dashed lines) with the experimental data. Main panel: am- 
plitude of the minimum of Av/v in function of temperature; 
inset: pressure variation for selected temperatures. Dashed 
lines show results assuming a gaussian distribution of Pc 

The pressure dependence of the dip is roughly repro- 
duced by this model but the agreement remains qualita- 
tive. In particular, as we approach Tc, the dip on the 
theoretical curves gets narrow much faster than the one 
observed experimentally. In terms of the scaling function 
/(x), the pressure profile is the most influenced by the 
behavior of / at intermediate values of x, where the scal- 
ing relation is not perfectly obeyed (for these simulations 
we took the average scaling function based on the curves 
reported in^). It is to note however, that independent of 
the scaling relations, such behavior is expected in our ap- 
proach as it is reminiscent of the pressure profiles of the 
conductivity. At any rate, the failure to reproduce the 
experimental curves with more accuracy cannot be con- 
sidered as a major shortcoming for this simplest model 
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where the complex effects of the strain on the multiband 
structure is parametrized with a simple scalar function 
D{V). 

In conclusion, we have studied the effects of the 
electron-lattice coupling in the vicinity of the critical 
point of the Mott transition line of Cr-doped V2O3. Crit- 
ical effects are only observed for the symmetry-invariant 
component of the strain tensor, which strongly supports 
the Ising scenario for this transition. The amplitude of 
the lattice softening is an order of magnitude smaller than 
estimations from a single-band Hubbard model in the 
DMFT framework. A simple model where the expression 
of the electron-lattice interaction is reduced to a simple 
scalar function D{V) gives a qualitative agreement and is 
consistent with the scaling relation for resistivity. How- 



ever a more realistic modeling would need to take into 
account the multi-orbital nature of this materia?^. 
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